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Motivated by the earlier experimental results and ab initio studies on the electronic structure 
of layered ruthenates (Sr 2 Ru 04 and Ca 2 Ru 04 ) we introduce and investigate the multiband d — p 
charge transfer model describing a single Ru 04 layer, similar to the charge transfer model for a 
single Cu02 plane including apical oxygen orbitals in high Tc cuprates. The present model takes 
into account nearest-neighbor anisotropic ruthenium-oxygen d — p and oxygen-oxygen p — p hopping 
elements, crystal-field splittings and spin-orbit coupling. The intraorbital Coulomb repulsion and 
Hund’s exchange are defined not only at ruthenium but also at oxygen ions. Our results demonstrate 
that the Ru 04 layer cannot be regarded to be a pure ruthenium t 2 g system. We examine a different 
scenario in which ruthenium Cg orbitals are partly occupied and highlight the significance of oxygen 
orbitals. We point out that the predictions of an idealized model based on ionic configuration (with 
no = 4-|-4x6 = 28 electrons per Ru 04 unit) do not agree with the experimental facts for Sr 2 Ru 04 
which support our finding that the electron number in the d — p states is significantly smaller. In 
fact, we find the electron occupation of d and p orbitals for a single Ru 04 unit n = 28 — x, being 
smaller by at least 1-1.5 electrons from that in the ionic model and corresponding to self-doping 
with X ~ 1.5. 

PACS numbers: 71.10.Fd, 71.70.Ej, 74.70.Pq, 75.10.Lp 


I. INTRODUCTION 

The description of the electronic states of transition 
metal oxides with partly filled d orbitals is not an easy 
task, and one usually looks for some simplifications. 
Models for manganites are complex as they involve partly 
filled t2g and Cg orbitals [T] . Systems with partly filled t2g 
orbitals and empty Cg orbitals, such as for titanates 
or vanadates [ 3 ] or with completely hlled t2g and partly 
filled Cg orbitals as in cuprates [ 4 ] or nickelates [ 5 ] are 
much easier to investigate. Such systems can be realized 
if the gap separating t2g from eg orbital states (induced 
by crystal-field effects) is sufficiently large. 

The transition-metal oxides with Ad ions are even more 
challenging as there electron correlations are somewhat 
weaker and simultaneously spin-orbit coupling plays an 
important role. Therefore one has to treat Ad electrons as 
both itinerant and strongly correlated, in the vicinity of a 
metal-insulator transition. The compound which belongs 
to this class and received a lot of attention is Sr2Ru04 as 
it became a candidate for a p-wave superconductor [^. 
Recent progress in photoemission technique made it pos¬ 
sible to investigate the many-body effects both in bulk 
and surface bands of Sr2Ru04 [ 7 ]. Spin-orbital entangled 
states have been seen recently in spin- and angle-resolved 
photoemission spectroscopy [ 5 ]. Such states arise in cor¬ 
related transition metal oxides either on superexchange 
bonds [ 9 ] or locally due to strong spin-orbit coupling at 
transition-metal sites [lO] . Indeed, these quantum effects 
play an important role in Mott insulators with Ad ions 
m and in these systems doped by 3 d elements [Hiig. 

In Sr2Ru04 the experiments support the earlier im¬ 
plementation of spin-orbit coupling within the local den¬ 


sity approximation with Coulomb interaction U treated 
in LDA-bC/ approach which show that both Coulomb U 
and the spin-orbit coupling are necessary for a correct 
description of the Fermi surface in Sr2Ru04 [ 13 ]. The 
electronic structure of Sr2Ru04 was extensively studied 
in the past m and it was established that the orbital 
physics plays here an important role | 16 j . In Ca2Ru04 
the bandwidth is smaller but the effects of spin-orbit cou¬ 
pling are even more pronounced m 

Due to a rather large crystal-field splitting between t2g 
to Cg levels ^ 3.5 eV in Ca2Ru04 and Sr2Ru04 [ 18 ], one 
might expect that these compounds are also purely t2g 
systems. This picture was also supported by the ear¬ 
lier studies of electronic structure in Ca2-a;Sra;Ru04 by 
photoemission m- Orbital polarization is then helpful 
to understand the x-ray absorption measurements and 
is of importance in describing the insulating state of 
Ca2-a;Sr2,Ru04 [SO]. A unique feature of Ca2-a;Sr2;Ru04 
is that slight changes in lattice parameters can induce 
drastic modifications of the character of their electronic 
ground states — Sr2Ru04 is metallic and superconduct¬ 
ing at low temperature [ 5 ], whereas Ca2Ru04 is distorted 
and undergoes a metal-insulator transition [ 21 ] . Interest 
in these materials is also motivated by very unusual invar 
effect reported recently [ 32 ] which suggests spin-orbital 
entanglement [S] in the ground state. 

Recently a simple tight-binding model was employed to 
investigate the superconductivity in Sr2Ru04 [ 33 ]. The 
simplest d — p model would just include three t2g or¬ 
bitals per Ru and three 2 p orbitals per O ion. However, 
there are serious doubts about the validity of this sim¬ 
plified physical picture. First, there are ah initio clus¬ 
ter-l-embedding computations by Kaplan and Soullard 
[ 31 ] who claim that the p-orbital charges on oxygens (in 
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Sr2Ru04) are not close to formal 6.0 as they follow from 
the ionic model, but are instead close to 5 . 0 , and that 
in addition Cg levels are partly occupied. At the same 
time the transfer of charge from any strontium ion (to 
ruthenium-oxide layer) is smaller than a formal value of 
2 electrons. Secondly, a similar system, namely C0O2 
layer was also believed to be pure t2g systems. However, 
it was shown | 25 j that Cg orbitals (in a doped system) can 
become very important, even in the absence of spin-orbit 
coupling. We show below that the total electron den¬ 
sity which follows from the idealized ionic model with 
formal electronic charges does not describe correctly the 
electronic structure of Sr2Ru04. In contrast, the d — p 
model with reduced electron density based on ab initio 
calculations | 24 j gives results which agree with experi¬ 
ment. 

To resolve the question about actual electron den¬ 
sity within eg orbitals in Sr2Ru04 we constructed a 
multiband charge-transfer model and performed unre¬ 
stricted Hartree-Fock computations on a finite RUO4 
cluster which contains 4 x 4 Ru ions and 4 x 4 x 4 oxy¬ 
gen ions — half of them located within the same plane 
as Ru ions, while the second half belonging to the elon¬ 
gated RuOe octahedra and surrounding the plane from 
above and below, being in out-of-plane (apical) positions. 
We imposed cyclic boundary conditions. The model in¬ 
volves (per a single RUO4 unit) five 4 d orbitals on Ru 
and 4 x 3 oxygen 2 p orbitals per unit cell occupied by: 
(i) no = 4 -I- 4 X 6 electrons, according to the formal and 
idealized ionic model; (ii) the electron number lower by 
at least one electron, i.e., n = 3 4 - 4 x 6 electrons or 
even smaller (we study below the case of n = 2 . 54 - 4 x 6 
electrons), as found in Ref. [ 24 ] . 

The paper is organized as follows. In Sec. jl^we intro¬ 
duce the multiband model which includes all Ad states at 
ruthenium ions and 2 p states at oxygen ions. The param¬ 
eters of the model are specified in Sec. |III[ The Hartree- 
Fock approximation for the Coulomb interactions is ex¬ 
plained in Sec. IV while in Sec. |^we present the results 
of numerical calculations and we introduce the concept 
of self-doping with respect to the electron densities in the 
ionic model. The paper is concluded with a short discus¬ 
sion and summary of the main results in Sec. VI In the 
Appendix we give the hopping elements d — p and p — p, 
respectively. 


II. MODEL HAMILTONIAN 


In this section we introduce the d — p charge- 
transfer Hamiltonian for RUO4 plane (such as realized 
in Sr2Ru04). It consists of several parts. 


T~L — Likin 4 - List 


H, 


diag 


Hi 


int ■ 


( 2 . 1 ) 


The different terms in Eq. ( 2 . 1 ) stand for the kinetic en¬ 
ergy (ilkin)) spin-orbit coupling {Hso), crystal-field split¬ 
tings which are diagonal in the {Ad, 2 p} orbital basis 


(LLdiag), and the intraatomic Coulomb interactions (Hint) 
— they all are explained below. 


A. Kinetic energy in hybridized d — p bands 
The kinetic part of the Hamiltonian is: 

Hkin = ^ ^ + H-C.^ , (2.2) 

where we employ a general notation, with ^ standing 
for the creation of an electron at site j in an orbital n 
with up and down spin, a =t, i- The model includes 
all Ad orbital states per Ru atom, v G {xy,yz,zx,x^ — 
y^, 3 z^ — r^}, and three 2 p orbitals per oxygen atom, u G 
{PxiPvtPz}- Alternatively, i.e., choosing a more intuitive 
notation, we can write dj ^ ^ for d orbitals, while p] £, 
for p orbitals. 

The matrix is assumed to be non-zero only for 

nearest neighbor ruthenium-oxygen d — p pairs, and for 
nearest neighbor oxygen-oxygen p — p pairs. The next 
nearest hopping elements, in particular direct ruthenium- 
ruthenium ones, and those between the p orbitals of 
neighboring apical oxygens are neglected. The nonzero 
ti g.j i, elements are listed in the Appendix. 


B. Spin-orbit coupling in layered ruthenates 


Formally, simplified spin-orbit part Hso of the Hamil¬ 
tonian Eq. ( 2 . 1 1 has a similar mathematical structure 
to the kinetic part Hkin PSH 5 D] . with elements 

restricted to single ruthenium sites. 




= E^io 


E 1(2.3) 

where the summation runs only over ruthenium sites and 
where we explicitly use d\ operators for the Ad or¬ 
bitals at Ru sites. The derivation of spin-orbit coupling 
starts from a single-site model. Using the {\i,pL,a){ 
basis one evaluates the full matrix of scalar products, 
{i, v, cr'|L.i • Si|z, fi, a), of angular momentum with spin 
Si operator (at site i). The individual single-site terms 
in are defined by on-site hopping-like elements, 

U v,(7'-,11,(7 (between different spin and o rbita l states), and 
one arrives at the matrix form in Eq. ( 2 . 31 . 


As we use here the basis of real Ad orbitals (and not the 
spherical harmonics) several spin-orbit elements turn out 
to be purely imaginary (hence the hermitian Hamiltonian 
is not real but complex). The elements of the matrix 
for a ruthenium site i are the following ones (for a similar 
result consult Fig. 6 in Ref. [ 27 ]), 
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where C is the spin-orbit coupling parameter and where 
the columns and rows are labeled in the following order: 

(xj/1) , {yz t) , {zx t), t), (3^^ - t) , 

(xy 4 .), {yz 4 .), {zx 4 .), - y^), (Sz^ - r^) ■ 

Note that the consequence of finite spin-orbit coupling 
C is nonconservation of the zth component of the total 
spin and therefore the obtained ground state wave func¬ 
tion is not a product of two Slater determinants for f- and 
4,-spin. In some cases the spin-orbit coupling Hgo can be 
treated as a minor perturbation and can be neglected 
when the average value of local angular momentum (at 
site i) is quenched to zero due to suitably strong crystal- 
field effects and low enough local symmetry (reduced by 
nearest neighbor atoms). Then, the spin-orbit Hamilto¬ 
nian can contribute to total energy only as second-order 
correction. Such a reasoning however allows one to make 
only qualitative predictions. 

An explicit treatment of spin-orbit coupling causes 
some difficulties. It is likely that the true ground states 
are not homogeneous in space, e.g. involving spin spirals 
or other micro-modulations, thus they may be considered 
intractable within Hartree-Fock computations for d — p 
clusters (they have too many order parameters to con¬ 
verge). One can only hope that these micro-modulations 
are of secondary importance. We assume this scenario 
and to have a tractable model we decided to use a sim¬ 
plified approach. Namely, we break the symmetry along 
natural quantization axis which is the z-th axis (perpen¬ 
dicular to Ru 04 layer). We emphasize that the aver¬ 
ages of local spin components aligned parallel to the (a, b) 
plane are assumed to be zero. To present this assumption 
in a more transparent way we can write down the formula 
for a local spin-flip, ^ ^ and a similar one 

for S~^. Thus the requirements that {S^^) = {S~^) = 0 
are equivalent to setting to zero the following order pa¬ 
rameters with opposite spins, {d] ^ ^d^ ^ _^) = 0. 

In Sec. |V] we report a study of charge space- 
homogeneous solutions. All occupation numbers, i.e., 
primary order parameters, are assumed to be the same 
for equivalent ruthenium ions and similar for oxygen 
ions. When studying the possibility of antiferromag¬ 
netism there are two sublattices — thus the number 
of order parameters doubles. Looking for charge space- 


r 


homogeneous ground states can be considered a simplifi¬ 
cation but from another point of view it can be treated 
as a consequence of strong long-range interionic electro¬ 
static interactions. These interactions are not explicitly 
included in the model (2.1) at present. We remark that 


long-range interionic interactions cannot be easily incor¬ 
porated into a typical d—p model but instead they can be 
accounted for by requiring that individual ionic charges 
are space-homogeneous in accordance with crystal sym¬ 
metry. 


C. Crystal-field splittings 

Let us now return to the d — p Hamiltonian ( |2.1[ ) . The 
next part of the model iLdiag is diagonal in the orbital 
basis and depends only on electron number operators. 
It takes into account the effects of crystal field and the 
difference of reference orbital energies (here we employ 
the electron notation), 

A = ed-ep, (2.5) 

between d and p orbitals, both for empty states without 
the Hartree terms which follow from iLint- Below we fix 
the reference energy Sd for d orbitals to zero, hence we 
use only Sp as a parameter and write: 

Hdiag = ^pPi,fj,,aPi,fi,cr 

i,fj.=x,y,z-,a 

+ E (2.6) 

i,^—xy,yz,...\(T 

Here the first sum is restricted to oxygen sites, while the 
second one runs over ruthenium sites. 

What concerns the value of Sp, it could be different for 
in-plane and for out-of-plane (apical) oxygens. Accord¬ 
ing to the earlier studies [SI1I31] this difference can be 
as large as 1.5 eV. But such a large value follows from 
a simplified procedure of fitting electronic bands to the 
LDA results. In the framework of the present d—p model 
one should expect much smaller difference, if any. The 
zero difference was assumed in computations performed 
in Ref. [33] and we also adopt this value. Here we remark 
that our test Hartree-Fock computations performed us¬ 
ing a big 1.5 eV splitting give large differences in charge 
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occupation between in-plane and apical oxygens, in dis¬ 
agreement with the results of the population analysis in 
Ref. [ 23 ]. This choice is indeed unrealistic as in addition 
one obtains only ferromagnetic ground states, see below. 

Let us return again to iLdiag- The vector containing 
the elements of can be expressed as 
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( 2 . 7 ) 


It includes the orbital splittings of Ad orbitals in the 
tetragonal crystal field. The constant Di serves as a 
crude estimate of the splitting between the orbital xy 
and the orbital doublet {yz,zx}, i.e., when taking into 
account only iLdiag while neglecting the remaining parts 
of the full Hamiltonian,. For Di < 0 (like in Ca2Ru04) 
in the presence of a tetragonal distortion of RuOe octahe- 
dra the xy orbital has a lower energy and is occupied (in 
the atomic configuration) by two electrons, while yz and 
zx are occupied by one electron each — then the energy 


gain is Hi, i.e., such an occupation pattern is more stable. 
For Hi > 0 (like in Sr2Ru04) the doublet {yz, zx} has a 
lower energy and is more stable. The parameters H2 and 
H3 follow from the estimates of the splitting between t2g 
and eg orbitals. They do not influence the ground state 
energy if only t2g orbitals are occupied, while there is a 
punishment (by the value of H2 or H3) for each electron 
occupying x‘^ — y^ or iz'^—r'^ level respectively. These pa¬ 
rameters {H2,H3} are much larger than Hi in the case 
of Sr2Ru04 and the difference between them is rather 
small. 

Jahn-Teller part can be neglected in the Hamiltonian 
( 2 . 1 ). The exception is the elongation of bonds between 
ruthenium and apical oxygens which could be considered 
as a frozen global (static) Qs Jahn-Teller distortion [T] 
but it is much simpler to include it by a proper renor¬ 
malization of the crystal-field splittings. Note that in the 
other compound Ca2Ru04 Jahn-Teller effects are large 
(see for example Ref. m) and as a result the symmetry 
of the corresponding RUO4 plane is significantly lowered. 


D. Local Coulomb interactions 

The last part of the multiband d — p Hamiltonian Hi^t 
stands for strong on-site Coulomb interactions. For the d 
orbitals at ruthenium sites it includes Hubbard intraor¬ 
bital repulsion Ud, Hund’s exchange and pair hopping 
also given by 






5 H 


fj,i> j 


( 2 . 8 ) 


Here is the tensor of Hund’s on-site interorbital exchange elements for d orbitals which can be expressed using 
Racah parameters B and C | 531 | 3 S] (see also Table I given by Horsch in Ref. [ 35 ]). Note that we sum twice over each 
pair {/ijz}) of orbitals in Eq. ( 2 . 8 ). Importance of local Coulomb interactions has been recognized in several model 
studies [ 37 | [ 38 |. In particular, strong correlations which originate from Hund’s coupling have been studied [ 39 ] and it 
was also suggested that this coupling supports the triplet superconductivity [ 30 ] . 

The anisotropy between different Hund’s exchange elements {may be neglected as long as one may limit oneself 
to the orbitals of the same symmetry, i.e., either to t2g or to Cg orbitals [ 9 ]. For convenience, we rewrite equation 
(2.8) as follows 


U^int ~ Ud 'y ] T o ^ ^ {Ud fiv) ^ ^ {Ud fiv) ' 




I,II 


Similarly, local Coulomb interactions at oxygen sites (for 2 p orbitals) are given by 

~ Up y ) T 2 ^ / ^Up ~ J^h) 'ni,p,rTni^u,a + ^ ) {Up — 2Jjj) 

i,li i,11:^1^,cr 

- yz + yz JHPlp,tPlp,lP^,l'APi,<^^ ^ 


( 2 . 9 ) 


(T't'Z,!/', —(T 




i,li^iy 


( 2 . 10 ) 
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where all off-diagonal elements are equal as they con¬ 
nect the orbitals of the same symmetry, i.e., = J^. 

Up to now this latter part neglected in the 

majority of studies, i.e., they assume Up = = 0. As 

a compensation some effective and appropriate modifica¬ 
tion of the charge-transfer energy has to be performed. 
Indeed, it has been shown for La 2 -a;Sr 2 ;Cu 04 cuprates 
employing constrained local density approach that the 
estimation of Up is very sensitive to the assumed value 
of A [?T] — the value of A = 2.0 eV yields Up ^ 8 eV, 
while A = 4.0 eV yields Up ^ A eV. One notes that there 
is a roughly linear dependence of Up on A to reproduce a 
constant value of the charge-transfer gap in Hartree-Fock 
[22]; at the same time the estimated value of Ud remains 
more or less constant. 

Unfortunately, not much is known about the real value 
of A (2.51 for different compounds. Estimates based on 
LDA results for Sr 2 Ru 04 gave the value A ~ 1.5 eV [52] . 
This value is much lower than A ^ 8 eV used in Ref. [53] 
(where also Up = 0 was assumed). All these estimations 
come from different fitting procedures employing numer¬ 
ous simplifying assumptions. One has to realize that A 
(2.5) is not an ab initio like value in the framework od 
d — p model but rather some effective value which better 
should be treated as a free parameter (the d — p model 
itself is an effective model, and definitely not an ab initio 
model). Thus in the following we will vary the value of A 
from 0 to 6.0 eV. We do not consider negative values of 
A as their consequence are greatly overcharged d-shells. 
We also remark that it has been realized nowadays that 
effective models with only d orbitals are insufficient and 
many papers treat electronic oxygen degrees of freedom 
explicitly. For cuprates this subject has a long history 
and realistic d — p models were studied for CuOa chains 
[421148| and for Cu02 planes [21 1221125] . The oxygen or¬ 
bitals are also of importance for ruthenates such as in 
Sr 2 Ru 04 [25] and for other correlated oxides 


III. SETTING THE HAMILTONIAN 
PARAMETERS 

A. Previous studies of the charge-transfer model 

The effective d — p model requires a choice of a num¬ 
ber of explicitly included parameters. In the extreme 
case when only Ad orbitals are used in a tight-binding 
(semiempirical) model [28] the parameters are very dif¬ 
ferent from the cases where Coulomb interactions are 
treated in the Hartree-Fock approximation. Here we 
adopt in-plane hopping elements (pda), {pdir), (jppa) and 
{ppir) used in [55]: —3.4, 1.53, 0.6, —0.15 (all in eV). The 
out-of-plane hoppings (involving the apical oxygens) were 
scaled using the formulae from the book by Harrison m- 
—2.6, 1.167, 0.559, —0.140 (all in eV). A similar value of 
(pdir) = 1.5 eV for the in-plane hopping was reported 
earlier by Oguchi [31] who applied the tight-binding for¬ 
mulae to fit the LDA electronic structure. Different es¬ 


timations for (pdir) are quite close to one another; they 
read as follows: (i) 1.0 eV in Ref. [5T]; (ii) ~ 1 eV as 
used by the group of Fujimori’s [30]; (hi) 0.85 eV in Ref. 
[55] . Let us note that for cuprates frequently accepted 
values for (pdir) are 0.75 eV in Ref. [2T] and 0.9 eV in 
Ref. [2]. 

The choice of the Coulomb elements is rather difficult 
due to their considerable screening in the oxides which is 
however less efficient in Ad systems [49] . For the intraor¬ 
bital Coulomb repulsion Ud at ruthenium sites the value 
of ~ 3 eV is most frequently used [SO] [351 [551 [53] . We fix 
here Ud = 3.1 eV following Ref. [55]. We remark that 
we cannot follow popular estimations made by Okamoto 
and Millis [52] and by Liebsch [55] as they apply to effec¬ 
tive models featuring only Ru sites renormalized by the 
hybridization with oxygen orbitals. 

Hund’s exchange elements are less screened than in¬ 
traorbital Coulomb elements and close to their atomic 
values (see for example Ref. [55]). For Hund’s exchange 
between two t 2 g electrons various estimates range 
from 0.5 eV up to 0.8 eV: (i) 0.5 eV in Refs. [331 [53]; (ii) 
0.6-0.8 eV in Ref. [57]; (hi) 0.7 eV in Ref. [55]; (iv) 0.8 eV 
in Ref. [TS]. We will use = 0.7 eV. Moreover, for the 
sake of fixing precisely Hund’s coupling tensor elements 
Jppj we use Table I from Ref. [36] and in addition we use 
an empiric formula C ~ AB for Racah parameters. With 
this Ansatz for a pure t 2 g system = 3B -|- C « 7i? 
and B = 0.1 eV [5]. This determines the ^ elements 
when Cg levels are not empty (using again the entries 
from Table I in Ref. [55]'). 

For the intraorbital Coulomb repulsion at oxygen sites 
Up (in ruthenates) again not much is known and it was 
neglected in several studies. In Ref. m this element is 
estimated to be Up = 4 — 6 eV. In cuprates the available 
data are more abundant: Up is 4.5 eV in Ref. [2], ^ 
4 eV in Refs. [m [55]; while several possibilities were 
also given (all in the range 3-8 eV) with 6 eV indicated 
by some experimental data [59]. We use below Up = 
4.4 eV. For Hund’s coupling at oxygen ions the values 
= 0.6 — 0.8 eV were suggested [57], while Grant and 
McMahan computations in cuprates yield = 0.8 eV 
[55] . Following these estimations, we use below = 0.8 
eV. Note that the corresponding tensor Jp has 

also the same entries for all off-diagonal elements. 

The spin-orbit coupling C on Ru sites is usually as¬ 
sumed to be in range from 0.10 to 0.17 eV [51 155113511331 
160] . Here we follow the most recent experiments which 
suggest that this coupling is in the middle of this range 
and take the value C = 0.13 eV [5]. 

To complete the set of the Hamiltonian parameters we 
have to provide estimates for the crystal-field splittings. 
The splitting between xy and {yz, zx} orbital levels is 
estimated as: (i) ~ 1 eV in Ref. [35]; (ii) ^ 0.3 eV in Refs. 
[35] [51]; (iii) 0.1 eV in Ref. [55]. We choose the value 0.1 
eV, i.e., we trust the reliable expertise presented in Ref. 
[55] . The splitting between t 2 g and Cg orbital levels is 3 
eV according to [T5] (in Ca 2 Ru 04 ) and up to 3.5 eV [30] . 
while the splitting of 0.8 eV between Cg orbital levels was 
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TABLE I. Parameters of the Hamiltonian (2.11 (all in eV) 
used for Hartree-Fock calculations. For the hopping integrals 
we adopt the values from |33| . Below we present only rep¬ 
resentative in-plane Slater integrals (pd-n) and (ppir). Out- 
of-plane integrals are obtained by applying Harrison scaling. 
Note that during computations we are setting the value of Ed 
to be zero as the reference energy. The value of A {—£p) is 
fixed to be the same for in-plane and for apical oxygens and 
is studied in the range A € [0.0, 6.0] eV. 


Ud 

jd 

Up 


C 

Dl 

D 2 

D 3 

(pdir) 

(ppiv) 

3.1 

0.7 

4.4 

0.8 

0.13 

0.10 

4.3 

3.5 

1.53 

-0.15 


assumed m- We accept these values in the parameter 
set employed in the present calculations (see Table I). 

Looking at the parameters, one important remark is 
proper. Namely the three different values: (i) the abso¬ 
lute value of parameter (pda) which is involved in hop¬ 
ping processes from t 2 g to Cg orbitals (see Table V in the 
Appendix); (ii) the value of Ud; and also (hi) the splitting 
between t 2 g and eg are all close to 3 eV. In other words 
— the splitting between t 2 g and eg does not seem to be 
large enough to justify the expectation that eg levels are 
almost empty. 


peat, once more, these results as they are really impor¬ 
tant for a proper understanding of the electronic struc¬ 
ture of Sr 2 Ru 04 : (i) the p orbital charge on oxygens (in 
Sr 2 Ru 04 ) is not formal 6.0 but is closer to 5.0 (oxygen s 
orbitals are also not fully occupied); (ii) the occupation 
on d orbitals is close to 6 but eg levels are partly occu¬ 
pied; (iii) charges on strontium ions are not formal 2 "’' 
but rather ~ 

Direct mapping of the ab initio results to the d — p 
model is not possible (as the d — p model neglects the 
valence s orbitals). However, it seems clear that the for¬ 
mal (idealized) ionic model with 6 electrons occupying p 
levels of each oxygen and 4 electrons occupying d levels 
of each Ru ion does not apply to the realistic Sr 2 Ru 04 . 
For the sake of convenience let as take a convention and 
introduce the self-doping x for a single RUO 4 unit with 
respect to the formal idealized model {x = 0 ) while for 
the real substance we shall consider finite self-doping val¬ 
ues such as 1.0, 1.25 and 1.5 (n = 28 — x, here we use 
such simple numbers so as the self-doping translates into 
integer electron number for the entire cluster). 

IV. THE UNRESTRICTED HARTREE-FOCK 
APPROXIMATION 


A. The self-consistent Hartree-Fock problem 


B. Motivation by earlier ab initio results 

Let us repeat that in the majority of the papers it is 
being assumed that eg levels are entirely empty. How¬ 
ever, in the present paper we make an attempt to deter¬ 
mine the electron densities in eg orbitals and to provide 
a realistic estimate of the charge on oxygens. We re¬ 
mind the reader that we are motivated by the ab initio 
computations performed on a small cluster-bembedding 
by Kaplan and Soullard |Mj. We take the liberty to re- 


We use the unrestricte d Ha rtree-Fock approximation 
to investigate the model (2.1) for Sr 2 Ru 04 . The tech¬ 
nical implementation is the same as described in Refs. 

1331 [i3]- Namely, the local Coulomb interaction 
Hamiltonian Hintra is replaced by Hartree-Fock mean- 
field terms. To give an example, the term according to a 
common interpretation of Wick, Bloch, and de Dominicis 
theorem the term Ud g be replaced (for 

Hartree-Fock computations) with one-electron operators 
and double counting correction terms. 


2,/J, 








■E 


i.M.t 


r)(4 








(4.1) 


Note that the terms with superconducting correlations 
are ignored in the above. Note also that standard mean- 
field decoupling usually ignores spin-flip terms (second 
line) [331132] • We remind that spin-fl ip te rms do appear 
in spin-orbit part of the Hamiltonian (2.1), see Eq. (2.3), 


and have to be included here on equal footing as the 
mean-field terms. 


The averages {dl ^ ^d^ ^ .^) and other similar ones can 
be treated as order parameters. At the beginning some 
initial values (a guess) have to be assigned to them. Dur¬ 
ing Hartree-Fock iterations the order parameters are re¬ 
calculated self-consistently until con vergence. When de¬ 
coupling all the terms in Hint ( 2 . 8 ) one finds that the 
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complete set of order parameters is as follows: 


where ^ ^ v. The off-diagonal elements (/r ^ v) are of 
crucial importance, particularly in the present case when 
finite spin-orbit coupling induces their finite values. A 
similar set of order parameters has to be considered for 
the p oxygen orbitals. 


B. Hartree-Fock calculations for Sr 2 Ru 04 


As we mentioned earlier, we are interested only in 
charge-homogeneous solutions, i.e., homogeneity con¬ 
cerns only primary order parameters (/r = v, i.e., charge 
occupations) but not off-diagonal (/r ^ v) order parame¬ 
ters. Thus primary order parameters are translationally 
invariant according to assumed lattice symmetry. In ad¬ 
dition, the symmetry is broken along the z-th axis and 
the order parameters {{d\ ^ ^ are fixed to be zero. 

Finally, the four-fold symmetry is imposed so that the 
occupation of Px and Py oxygen orbitals is the same, and 
also of yz and zx orbitals for Ru ions. Altogether we 
have got 7 independent primary order parameters (per 
Ru 04 unit) when looking for paramagnetic ground state 
and 15 for ferromagnetic (or antiferromagnetic) ground 
states. These numbers are large enough to expect trou¬ 
bles with the Hartree-Fock convergence and indeed this 
is the case. The regular convergence is found for a very 
limited set of the Hamiltonian parameters (for example 
for situations when the oxygen occupations are very close 
to 6 as was the case in Ref. [33]). A typical situation for 
our computations is that Hartree-Fock iterations do not 
converge but oscillate (in a two-cycle) instead. 

The standard remedy for poor convergence is the so- 
called dumping technique, but here it failed, unfortu¬ 
nately. We had to resort to quantum chemistry tech¬ 
nique called level shifting [^. It is based on replacing 
the Hartree-Fock Hamiltonian by a different Hamiltonian 
— the one with the identical eigenvectors (one parti¬ 
cle eigenfunctions) and with identical occupied eigenen- 
ergies. The original eigenenergies of virtual states are 
all uniformly shifted upwards by a fixed constant value. 
Thus if we apply the shift say by 1 eV, then the HOMO- 
LUMO gap (the gap between highest occupied and lowest 
unoccupied eigenstate) we obtain will be artificially en¬ 
larged exactly by I eV. 

When applying virtual level shifting technique we can 
obtain valuable information. First case is when the 
HOMO-LUMO splitting (after correcting for the shift) 
is negative (for a few different shifts and a few differ¬ 
ent starting conditions). Then the single-determinant 
Hartree-Fock ground state we obtain is probably not 
correct and single-determinant description of the ground 



FIG. 1. A fragment of the studied Cu 04 cluster with four Ru 
ions shown by dark (red) big spheres and surrounding them 
O ions shown by smaller gray (dark blue) spheres. The cal¬ 
culations in Sec. V were performed for a 4 x 4 Ru 04 periodic 
supercell. 


state is not possible at all. Multi-configuration Hartree- 
Fock is required in such a situation instead (and let us 
remind that the multi-configuration Hartree-Fock still 
did not mature enough to be a standard working tool in 
solid state physics). Second case is when HOMO-LUMO 
gap we obtain is zero (or very close to zero). Then the 
ground state identification is questionable. However, the 
probability that such an identification is correct and that 
the ground state is conducting can be substantial. The 
probability of correct identification can be further en¬ 
hanced when performing numerous extra computations: 
if we obtain the same identification for different shifts 
and different starting Hartree-Fock conditions, the result 
is accepted. Finally for the cases with positive HOMO- 
LUMO gap we usually had no problems. 

The computations were performed for a periodic 4x4 
Ru 04 supercell with its fragment shown in Fig. a (for 
each particular set of Hamiltonian parameters). They 
were repeated many times for different starting condi¬ 
tions (different starting charge occupations) and shifts. 
The state with the lowest Hartree-Fock energy was then 
identified as the ground state. Numerous runs were nec¬ 
essary as the Hartree-Fock convergence provides many 
different metastable states (metastable, i.e., only local 
but not a global minima of the energy). Typical num¬ 
ber of runs should be large enough — in some situations 
more than a hundred. For this very reason the detailed 
investigation of the phase diagram is too expensive (even 
for such a small cluster size as we use). Still we per¬ 
formed as many computations for as many different sets 
of Hamiltonian parameters so as to be sure about the 
general trends occurring on the phase diagram which are 
presented in Sec. V. 


V. NUMERICAL RESULTS 


Numerical studies of the multiband d — p model (13 
require not only the parameters which were fixed in 


Sec. Ill but also an assumption concerning the to¬ 
tal electron number per unit cell. We consider below 
two different scenarios: (i) the formal ionic model with 
710 = 4-1-4x6 = 28 electrons per RUO 4 unit, and (ii) the 
model with a smaller total number of n = 28 — a: elec¬ 
trons, where we investigated a few representative values 
of self-doping x = 1.0,1.25,1.5. Thereby we concentrate 
on the most important results obtained for a realistic 
value of Coulomb interaction within oxygen 2p orbitals. 
Up = 4.4 eV, within the framework of these two different 
scenarios. As we shall show below, these two situations 
require quite different values of Coulomb parameter Up 
at oxygen 2p orbitals. 


A. First scenario: Formal ionic model 

Taking the ionic model as a starting configuration for 
the Hartree-Fock iterations, we assume that each oxygen 
ion has 6 electrons within 2p orbitals and each Ru"*^'*' 
ion has 4 electrons within Ad orbitals. The RUO 4 unit 
has a negative charge Q = Ae which is compensated by 
two Sr 2 + ions considered only as electron donors to the 
plane of RUO 4 units. This charge distribution (assumed 
on start of Hartree-Fock iterations) is however unstable 
and the electrons quickly redistribute along the iteration 
process due to d—p hybridization. The final charge distri¬ 
bution is shown in Fig. [^and in the upper part of Table 



FIG. 2. Electron densities in the multiband model ( |2.1[ | for 
increasing A (A = —Sp) obtained in the formal ionic model 
(with no = 28 electrons per RUO 4 unit): at Ru ions (solid 
lines) and at O ions (dashed lines). Data points show electron 
densities within t 2 g orbitals {nt 2 g, circles), all 4ci orbitals (n 4 d, 
triangles), 2 p orbitals at in-plane oxygens (np\\, -f), and 2 p 
orbitals at apical oxygens {np±, x). The obtained ground 
state is ferromagnetic. 


TABLE II. The ground state obtained in the formal ionic 
model with no = 28 {x = 0 ) electrons and in the realistic 
model with self-doping (x > 0): HOMO-LUMO gaps G (eV), 
the total magnetic moment mtot per RUO 4 unit, and the mag¬ 
netic moment at Ru ions ttiru obtained for several values of 
A (eV). Note that when the total magnetization mtot is large, 
it is mainly due to the magnetization at oxygen ions. 



energies (eV) 

magnetizations 

X 

A 

G 

mtot 

mRu 

0.00 

0.0 

0.10 

1.00 

0.34 


1.0 

0.29 

0.75 

0.32 


2.0 

0.39 

0.75 

0.36 


3.0 

0.43 

0.75 

0.40 


4.0 

0.46 

0.75 

0.45 

1.00 

0.0 

0.27 

1.25 

0.30 


1.0 

0.10 

0.24 

0.09 


2.0 

0.30 

0.24 

0.11 


3.0 

0.33 

0.24 

0.12 


4.0 

0.35 

0.24 

0.14 

1.25 

0.0 

~ 0 

1.12 

0.26 


1.0 

~ 0 

0.12 

0.11 


2.0 

0.04 

0.11 

0.05 


3.0 

0.04 

0.11 

0.06 


4.0 

0.04 

0.11 

0.06 

1.50 

0.0 

0.22 

0.20 

0.18 


1.0 

0.21 

0.00 

0.00 


2.0 

0.26 

0.00 

0.00 


3.0 

0.38 

0.00 

0.00 


4.0 

0.51 

0.00 

0.00 


II. As expected, the total ( 714 ^) and partial {nt 2 g) elec¬ 
tron densities at Ru ions increase with increasing value 
of Ep (i.e., decrease with increasing value of A) which 
follows from electron transfer from O to Ru ions. 


The ground state (in the entire range of the investi¬ 
gated values of Ep) is ferromagnetic and insulating, which 
is surprising and does not agree with ab initio calcula¬ 
tions m- It could be argued that Hund’s exchange is 
strong enough to polarize the Ad electrons if t 2 g orbitals 
are well away from half filling {nt 2 g > 4.7 for the consid¬ 
ered range of Ep), in spite of large d — p hybridization. 
The instability towards ferromagnetism competes here 
with d — p hybridization and therefore the magnetic mo¬ 
ment per Ru 04 unit is small. The nonmagnetic ground 
state is metastable, and has a higher energy by about 0.8 
eV (per RUO 4 unit). Also, the occupation patterns for 
the obtained ground states do not agree with the ab initio 
data [24]. It is well known that Sr 2 Ru 04 has a paramag¬ 
netic and metallic ground state. Thus, the model (2.1) 
for the adopted parameter values and within the first 
scenario is clearly not realistic enough for Sr 2 Ru 04 . 
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B. Second scenario: Realistic self-doping model 

In the second scenario we follow the results presented 
in [23] and we assume a reduced total number of elec¬ 
trons per Ru 04 unit. Taking the total electron number 
n = 28 — X with x > 0 this corresponds to finite hole 
self-doping, and we study here x = 1.00,1.25,1.50. For 
the missing electrons the smaller (than in the ideal-ionic- 
model) transfer of valence 4s electrons from Sr sites is 
mainly responsible (5s valence electrons on Ru, neglected 
in the d — p model, have also some minor influence). The 
corresponding electron densities obtained for x = 1.50 
are shown in Fig. [^and in Table II. 

The difference between the two density distributions 
shown in Figs. 2 and 3 is mainly visible in the electron 
densities at Ru ions. At hnite self-doping of x = 1.5 
the total electron density within 4d orbitals is close to 
n-id = 5 for a value A ~ 5.0 eV, while it is close to 
nid = 5.5 for the same value of A in the ionic model. 
The electron density at in-plane oxygens is also somewhat 
reduced in the former case. These orbitals are influenced 
stronger by the self-doping as they are hybridized with 
the t 2 g orbitals at the central Ru ion in each RUO 4 unit, 
and provide also bonding between these units. 

Our calculations demonstrate that the metallic non¬ 
magnetic state observed in Sr 2 Ru 04 is realized when the 
density of electrons within Ru t 2 g orbitals is not close 
to half filling, i.e., nt 2 g < 4.5. Such a nomagnetic state 
is found in the present realistic model with self-doping 
when of 2p energy A > 1.0 eV. On the contrary, when 
X G [1.25,1.50] and A G [0.0,1.0] eV, the ground state is 



FIG. 3. Total electron densities in the realistic d — p multi¬ 
band model based on ab initio calculations for increasing A 
(A = —£p) for Ru(t 29 ) orbitals (solid line, circles), Ru(4d) or¬ 
bitals (solid line, triangles) and 0 ( 2 p) orbitals (dashed lines). 
There are (4x6-1- 2.5) electrons per a single RUO 4 unit, cor¬ 
responding to self-doping by x = 1.5 electrons. The ground 
state is nonmagnetic for A > 1.0 eV, and weakly ferromag¬ 
netic for A < 1.0 eV. The change of magnetic order is respon¬ 
sible for the change of the slope of the lines. 


TABLE III. Hartree-Fock energy Ahf, the HOMO-LUMO 
gap G, electron densities within t 2 g orbitals (nt 2 g), for all 4d 
orbitals ( 7144 ), oxygen 2p orbitals at in-plane oxygens (np||), 
and oxygen 2p orbitals at oxygens in apical positions (npj_), 
and total magnetic moments mtot (per RUO 4 unit), for a 
few selected ground states with the total electron number 
n = 28 — X per RUO 4 unit. Parameter: A = 1.0 eV. 

energies (eV) electron densities 


X 

Fhf 

G 

nt 2 g 

nid 

Up 11 

?lp L 

mtot 

1.00 

111.049 

0.24 

4.78 

6.19 

4.95 

5.45 

0.10 

1.25 

107.838 

0.12 

4.85 

6.21 

4.95 

5.32 

0.00 

1.50 

104.663 

0.00 

4.77 

6.13 

4.89 

5.30 

0.21 


ferromagnetic. Metallic ferromagnetism is here possible 
due to large values of the Stoner parameter I ~ Ud + ‘2.Jfj 
for partly filled t 2 g orbitals [52| which is farther enhanced 
by partly occupied Cg orbitals. Effectively I is enhanced 
here when the density of 4d electrons is increased and Cg 
orbitals are also partly filled as we have seen in the ionic 
model. 

The most interesting data obtained in our Hartree- 
Fock calculations for the realistic model are presented in 
Table III. Here we identify the region of phase diagram 
where both the HOMO-LUMO gap G is small or vanishes 
and the gronnd state is close to a transition from a non¬ 
magnetic to ferromagnetic one. In this regime one hnds 
large electron density within Cg orbitals accompanied by 
rather strong reduction of electron density at the oxy¬ 
gen ions in Ru 02 planes. The density at these ions np\\ 
varies from 4.95 to 4.89 when x G [1.0,1.5], i.e., each oxy¬ 
gen ion contains one hole and is rather close to the 
ionic state. These results of our computations agree well 
enough with the results of Ref. [24]. Indeed, the charge 
on oxygen ions, in particular the ones lying within RUO 2 
planes, is close to 5 and not to formal 6 electrons per ion. 
This demonstrates the metallic character of the electronic 
structure in Sr 2 Ru 04 . We emphasize that the occnpa- 
tions which follow from the present model are close to 
those reported in the ab initio investigation |24] — some 
representative examples of the occupations obtained in 
the Hartree-Fock calculations are shown in Table HI (see 
also Fig. 3). 

C. Importance of Coulomb interactions at oxygen 

ions 

We have verified that actual electron densities within 
2p orbitals are rather sensitive to the used Coulomb inter¬ 
action parameters at oxygen ions. To obtain nnphysical 
density of ripn ~ 6 at in-plane oxygen ions (as in the for¬ 
mal ionic model) one must require that Up = 0 and/or 
£p must be very large negative, i.e., the charge-transfer 
gap A has to be very large. This is confirmed by all test 
computations performed varying the values of Up and Sp. 

Using the set of Hamiltonian parameters from Table 
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TABLE IV. Hartree-Fock results for ferromagnetic ground 
states obtained using three different miltiband models in ab¬ 
sence of electron interactions at oxygen ions {Up = = 0): 

total electron densities within t 2 g orbitals {nt 2 g), for all Id or¬ 
bitals {n 4 d), oxygen 2 p orbitals at in-plane oxygens (u.p||), and 
oxygen 2 p orbitals at oxygens in apical positions {np±), to¬ 
tal magnetization nitot per RUO 4 unit, HOMO-LUMO gap G 
(eV), and spin-orbit contribution per single Ru (L/io ) (eV). In 
the third segment all HOMO-LUMO gaps G are zero (within 
1 meV accuracy), thus the reliable identification of the ferro¬ 
magnetic ground state is not possible (but probable). 




energies (eV) 


electron densities 

model 

X 

A 

G 

{X:j) 

nt2g 

Ueg 

U-pII 

?T-p L 

mtot 

I, all 4d 

0.0 

0.0 

0.39 

-0.076 

4.19 

0.51 

5.75 

5.91 

1.0 



3.0 

0.34 

-0.083 

4.12 

0.40 

5.81 

5.93 

1.0 



6.0 

0.16 

-0.088 

4.08 

0.31 

5.85 

5.95 

1.0 

I, only t 2 g 

0.0 

0.0 

0.37 

-0.043 

4.24 

- 

5.91 

5.98 

1.0 



3.0 

0.35 

-0.049 

4.14 

- 

5.94 

5.99 

1.0 



6.0 

0.20 

-0.052 

4.09 

- 

5.96 

5.99 

1.0 

realistic 

1.5 

0.0 

~ 0 

-0.026 

2.90 

0.62 

5.63 

5.86 

1.25 



3.0 

~ 0 

-0.025 

2.76 

0.50 

5.73 

5.91 

1.25 



6.0 

~ 0 

-0.025 

2.67 

0.37 

5.80 

5.93 

1.25 


1 but setting Up = = 0 (while keeping other param¬ 

eters unchanged) we performed additional Hartree-Fock 
computations to investigate importance of Coulomb re¬ 
pulsion at oxygen ions. These calculations gave very dif¬ 
ferent electron density distributions than those obtained 
before for the same values of A and total electron den¬ 
sity n, but with finite Up and J^. As shown in Table 
IV, one finds large electron densities at in-plane oxy¬ 
gens, ripii > 5.75, and only ferromagnetic ground states 
in the entire range of Sp G [—6,0]. It may be consid¬ 
ered quite unexpected that the 2 p oxygen orbitals are 
almost completely filled then even at Sp — 0. The al¬ 
most ionic state is here a consequence of Coulomb 
repulsion at Ru ions which blocks electron redistribution 
due to hybridization. It is also surprising that the same 
Hartree-Fock energy is obtained for ferromagnetic and 
for antiferromagnetic ground states. Note that for the 
pure t 2 g model (when setting II 2 3> 1 and Da S> 1 ) and 
for large A = 6.0 eV (when t 2 g occupation number is 
equal to formal nt 2 g = 4) the nonmagnetic ground state 
has a higher energy by ^ 0.8 eV than the ferromagnetic 
one. 


VI. DISCUSSION AND SUMMARY 

Altogether the results of the presented Hartree-Fock 
computations are too complex to be fully conclusive, 
but nevertheless this study uncovers several important 
facts concerning the modeling of ruthenium oxides by the 
multiband charge-transfer model. First of all, the d — p 
model with a minimal basis set consisting of {Ad, 2p} 


orbitals is a useful tool for investigating the electronic 
structure when the effective electron density within the 
considered basis set is established in agreement with the 
experimentally observed ground states. We have found 
that a significant electron charge is transferred beyond 
the d — p orbitals and thus the effective electron density 
within the RUO 4 units has to be reduced to n = 28 — a:, 
with X G [1.0,1.5]. This effect is similar to the reduction 
of electron density in d — p orbitals in cuprates, where 4s 
orbitals at Cu ions are also partly occupied [63] . Another 
reason responsible for this appreciably reduced electron 
density could be a partial charge transfer from oxygen or¬ 
bital to the charged Sr ions in Sr 2 Ru 04 , suggesting that 
the ionic picture with Sr^"*" ions transferring 2 electrons 
to the RUO 4 subsystem is oversimplified. 

We have shown that Coulomb interaction effects at 
oxygen ions are very important and have to be included 
in a realistic description of these materials. Only then 
the hybridization effects are strong enough and are able 
to overcome Hund’s exchange at ruthenium ions which 
otherwise induces metallic ferromagnetic state, contrary 
to the experimental observations. We note however that 
ferromagnetic instability was observed in Ca 2 -a;Sr 2 ;Ru 04 
systems where antiferromagnetic interactions are also 
possible |64j . This is reminiscent of the situation encoun¬ 
tered in the Cai_a;Sra;Ru 03 perovskites, with CaRuOa 
found to be on the verge of a ferromagnetic instability 
[65] . Therefore, we suggest that further research in a 
model including lattice distortions is required to estab¬ 
lish the range of stability of ferromagnetism in the lay¬ 
ered Ca 2 -a;Sra;Ru 04 systems, being to some extent also 
expected from the present results. 

An extension of this model could be used for a simi¬ 
lar modeling of the electronic structure of Ca 2 Ru 04 , but 
this would also require Jahn-Teller coupling to the lattice 
to describe correctly the lattice distortions which accom¬ 
pany the insulating state. It remains a challenge for the 
theory to establish whether the electron density within 
the multiband d—p model would not be increased by such 
an insulating state and we suggest that the self-doping 
effect described here would be concentration dependent 
in the Ca 2 -a;Sr 3 ;Ru 04 compounds. 

In summary, the most important consequence of both 
d—p hybridization and spin-orbit coupling is the electron 
transfer from t 2 g to Cg orbitals as our calculations demon¬ 
strate that eg orbitals at Ru ions are partly occupied in 
the realistic parameter regime. Thus, Cg orbitals have to 
be included in any realistic model for ruthenium oxides. 
This invalidates the paradigm that ruthenium oxides are 
pure t 2 g systems. A second very important effect is a 
significant reduction of the electron density within oxy¬ 
gen orbitals from the values obtained in the ionic model, 
which effectively corresponds to one hole per oxygen ion 
within Ru 02 planes. Finally, only when the above self¬ 
doping effect is included, the nonmagnetic metallic state 
of Sr 2 Ru 04 may be correctly described. 
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TABLE V. The non-zero ruthenium-oxygen hopping elements 
in Ru 04 plane as obtained using Slater-Koster rules [501166j . 
{pdn) and {j>da) are the appropriate Slater-Koster interatomic 
integrals [501 166] . We assume that oxygen ions belonging 
to the RuOe octahedron are reached by in-plane vectors 
±ai = ±a(l,0,0), ±a 2 = ±a(0,1,0), and the apical oxygen 
positions by out-of-plane vectors, ±a 3 = ±6(0,0,1). We use 
standard notation, with R = Ri — Rj being the vector for a 
nearest neighbor bond between ruthenium at site i and oxy¬ 
gen at site j, while the coordinates (l,m,n) of R/i? which 
stand for the direction cosines of the hopping. 


R 

{1, m, n) 

V 

M 


±ai 

(±1,0,0) 

X 

2 2 

X -y 

K^)ipdf^) 


(±1,0,0) 

X 

3z^ - 

-l{\){pda) 


(±1,0,0) 

V 

xy 

l{pdn) 


(±1,0,0) 

z 

zx 

l{pd-K) 

±a2 

(0,±1,0) 

y 

2 2 

X -y 

-m{^){pda) 


(0,±1,0) 

y 

3z^ - 

-m(|)(pdo-) 


(0,±1,0) 

X 

xy 

m{pdn) 


(0,±1,0) 

z 

yz 

m{pdn) 

±a3 

(0,0, ±1) 

z 

3z^ - 

n{pda) 


(0,0, ±1) 

X 

xz 

n{pd'K) 


(0,0, ±1) 

y 

yz 

nijpd-R) 
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APPENDIX 


The non-zero ruthenium-oxygen d — p and oxygen- 
oxygen p — p hopping elements obtained by using Slater- 
Koster rules [55] for the lattice constant a = 1 are pre¬ 
sented in Tables V and VI, respectively. 
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